function dydt = ED_fixedbed(t, y)
global NE NP ffd cf cf1 cf2 dpp AZ BZ lgt

C = y(1:lgt); %Salz
C1 = y(lgt+1:2*lgt); %Protein 1
C2 = y(2*lgt+1:3*lgt); %Protein 2
[dCdt, dC1dt, dC2dt] = ED_Bandfun1(t, ffd, C, C1, C2, NE, NP, cf, cf1, cf2, dpp, AZ, BZ);

dydt = [dCdt; dC1dt; dC2dt];